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ABSTRACT 


Unlike stochastic foams, 3D printed lattice structures can be created with a tailored microstructure to achieve the 
desired global mechanical properties. In this study, quasi-static energy absorption of six polymeric lattice struc- 
tures of different unit cell topologies and deformation behaviours but with identical relative density are investi- 
gated. For the first time, a recently developed powerful 3D printing technology, namely Multi Jet Fusion (MJF), is 
utilized to produce high quality and performance lattices. Nonlinear finite element analysis (FEA) is used to an- 
alyse their compressive response and energy-absorbing characteristics. To validate the FEA results, experimental 
compression tests are performed on two types of lattices. The results of this investigation showed that the global 
energy absorption of lattices can be improved through controlling and manipulating their micro-topology. It is 
found that the lattice structures with bending-dominated deformation have low stiffness and strength but pro- 
vide good energy absorption capability. The stretch and buckling dominated structures are shown to be stiffer 
and stronger but have low energy absorption performance. The findings of this research identify high performing 


Energy absorption unit cell geometries that can be used for energy absorption application of lattice structures in the development of 
advanced 3D printed structures that are superior to stochastic foams. 

© 2018 Elsevier Ltd. All rights reserved. 

1. Introduction undergo a large compressive strain at nearly a constant stress level, 


One of the major uses of cellular materials is in energy absorption 
applications, such as packaging and protective devices. They can 
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thus absorbing a large amount of energy without producing a high 
stress on the protected object. This makes them excellent for such appli- 
cations because the aim of using energy-absorbing materials is not just 
to absorb a certain amount of energy, but also to keep the transmitted 
force to the protected object below the threshold that can cause damage 
or injury [1]. Foams are the most well-known type of cellular materials 
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Table 1 


Unit cell and overall specimen design and dimensions for the examined 3D structures of 15% relative density. 


Structure type Lattice designation Unit cell design Cell design A2 x 2 x 2 cells structure showing Overall cubic spacemen 
parameters the way the cells joined together 
Dorl d CAD model Side length 
(mm) (mm) (mm) 

1 Circular D=10 1.488 51.488 

2 Octagonal l= 4.14 1.626 51.626 

3 Strengthened Octagonal l= 414 1.214 51.214 

4 Kelvin l= 3.54 1.496 51.496 

5 Rhombicuboctahedron (RO) 1=414 1.134 51.134 

6 Cubic 1=10 2.246 52.246 


that are widely used. Polymeric and metallic foams are tailored for a 
specific energy absorption application by choosing the right material 
and density. However, the drawback of foams is they are stochastic 
and have a random microstructure [2]. As the structure of these mate- 
rials at the micro level play a big role in their global behaviour and prop- 
erties [3,4], researchers have tried to find better alternatives for foams. 
Lattice materials (structures) fabricated by 3D printing can overcome 
this limitation to a great extent. By exploiting this advantage offered 
by 3D printing technology, lattice materials can be made with uniform 
and ordered microstructure, their unit cells can be manipulated and op- 
timized to attain the desired mechanical properties for a particular ap- 
plication [2]. 

Lattice structures are commonly constructed by duplicating three- 
dimensional mesoscale unit cells. They are multi-functional materials 
that can offer a range of desirable properties. They offer the potential 
to easily engineer specific meso-scale properties (at cell level) to pro- 
duce desirable macro-scale material properties for a wide variety of en- 
gineering applications [5]. 

Several researchers have investigated the influence of unit cell shape 
on compressive strength of lattice structures especially for bio-based 
tissue engineering (TE) scaffolds [6-10]. Others have studied the effect 
of the unit cell on the global deformation and crushing behaviour of cel- 
lular structures. For instance, Li et al. [11] showed that the unit cell 
shape design is an effective way to control the mechanical properties 
of the reticulate meshes such as the elastic modulus, compressive 
strength and deformation behaviour. Ashby [3] considered the distinc- 
tion between stretch-dominated and bending-dominated structures 
the most important concept in analysing the mechanical behaviour of 
lattice materials. For a given density, the first is extremely stiff and 
strong, while the second is compliant and not strong, but absorbs en- 
ergy well during compression. Ahmadi et al. [12] experimentally 


showed that the mechanical properties, compressive behaviour and 
failure mechanism of 3D printed porous titanium alloy biomaterials 
were highly dependent upon the type and dimensions of the building 
unit cells. Craddock et al. [13] indicated that despite the current limita- 
tion of additive manufacturing techniques they can be used for 
manufacturing tailored lattice structures for impact absorption. 
Vesenjak et al. [14] experimentally studied the influence of cell shape 
on the compressive behaviour and energy absorption of 2D polymeric 
cellular structures with quadratic and circular cells under quasi-static 
and dynamic uniaxial compressive loading. Hammetter et al. [15] inves- 
tigated the effect of strut slenderness, strut inclination angle, and num- 
ber of repeat lattice layers on the mechanics of compressive 
deformation including the specific strength and energy absorption of 
pyramidal lattice structures using analytical and numerical methods. 
Ozdemir et al. [5,16] experimentally and numerically studied the energy 
absorption of additively manufactured metallic lattice structures of 
cubic, diamond and re-intrant lattice structures under quasi-static and 
dynamic loading. Ullah et al. [17] used single unit cell approach to pre- 
dict the energy absorption capacity of Kagome and atomic truss config- 
urations built through selective laser melting via experimental and 
numerical studies. Gautam et al. [18] studied the influence of build ori- 
entation, truss radius and surface roughness on stiffness, strength and 
energy absorption of ABS Kagome truss single unit cells built via fused 
deposition modeling (FDM). Ha et al. [19] designed and studied an en- 
ergy absorbing lattice constructed from multiple tetra-beam-plate unit 
cells with negative stiffness behaviour. The lattice was fabricated by se- 
lective laser sintering (SLS) approach and analysed both numerically 
and experimentally. Kaur et al. [20] investigated the deformation of 
two types of 3D printed stretch-dominated micro-lattice structures 
made of different polymeric materials using compression testing and 
FEA simulation methods. Choy et al. [21] investigated the mechanical 
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Fig. 1. Schematic of building process in Multi Jet Fusion (MJF) system. 


properties of density graded metallic lattice structures of cubic and hex- 
agonal unit cells fabricated by selective laser melting (SLM). Al-Saedi 
et al. [22] studied the compressive properties and energy absorption ca- 
pability of functionally graded aluminum alloy lattice fabricated by SLM. 
Harris et al. [23] experimentally investigated the impact response of ad- 
ditively manufactured stainless steel hybrid lattice materials. Winter 
et al. [24] numerically studied the dynamic energy absorption behav- 
iour of a series of two-dimensional stainless steel cellular structures of 
different architectures. Mohsenizadeh et al. [25] experimentally 
showed that it is possible to design and additively manufacture light- 
weight polymeric metamaterials for energy absorption capable of re- 
covering their shape after high deformation. 

While most of the published works have investigated the behaviour 
of specific types of individual lattice structures with specific unit cell ge- 
ometry, very few studies have been made to compare the mechanical 
response and energy absorption characteristic of 3D printed lattices of 
different types of unit cells made of identical relative density. Re- 
searchers have also been able to design and 3D print energy absorbing 
lattice structures, but still more detailed studies are required to identify 
the idealised geometry of the unit cell, which will generate a lattice 


capable of absorbing maximum energy at the lowest transmitted stress 
level. 

This work examines the effect of a variety of 3D unit cell topologies 
on the energy absorption behaviour and compressive response of the 
respective 3D printed polymeric lattice structures of uniform relative 
density. The aim is the development of a high energy-absorbing lattice 
with the closest behaviour to the ideal energy absorber (the highest 
specific energy absorption at lowest stress level). Full-scale nonlinear fi- 
nite element analysis (FEA) method was used to simulate their com- 
pressive behaviour under uniaxial quasi-static loading. Experimental 
compressive tests are performed to validate the FEA results with lattice 
structures printed on a newly developed powerful 3D system, namely 
Multi Jet Fusion (MJF), in a new brand of polyamide 12 material namely 
HP PA12A available on the system. A systematic approach is presented 
to characterise and investigate the energy absorption properties of the 
3D lattice structures and a clear criterion is set to identify the best struc- 
tures. The global stress-strain behaviour and the energy absorption per- 
formance of the structures are related to the deformation behaviour of 
their building unit cells. The optimum lattice for energy absorption ap- 
plications is identified. 
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Type 5- Rhombicuboctehedron (RO) 
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Fig. 2. 3D printed prototypes of the lattice structures on HP Jet Fusion 4200 3D printing machine in polyamide 12 (HP PA12). 


2. Design and material 
2.1. Lattices design 


The relative density is the single most important feature of cellular 
solids [1]. It is the ratio of the apparent density of a cellular material, 
p*, to that of the parent solid material from which the cell edges are 
made, ps. It is also equal to the volume fraction of solid in the cellular 
material. In this study, to investigate the compressive response and en- 
ergy absorption performance of ordered lattice structures, this parame- 
ter was kept constant. 
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Fig. 3. Stress-strain curves obtained from standard tensile testing of HP PA12 coupons 
(specimen Type I, ASTM D638-14). 


Six types of cellular structures were designed from different unit 
cells. The unit cell shapes were chosen in a way that provides different 
deformation mechanisms of the cell edges, such as bending, stretching 
and buckling. Each unit cell had centre line dimensions of 10 x 10 
x 10 mm? and was geometrically symmetric around the three perpen- 
dicular planes intersected at its central point. Each lattice structure 
consisted of 5 x 5 x 5 cells constructed by duplicating unit cells in the 
three directions. The six structures were designed to have a constant 
relative density of 15%. Table 1 shows the computer-aided design 
(CAD) models of the six building unit cells, the way they were joined to- 
gether, the 3D lattice specimens constructed from each unit cell, and 
their dimensions. It can be observed from Table 1 that the unit cells 
that have more edge members have their struts' diameter thinner in 
order to keep the relative density constant. 

The unit cell of type 1 (Table 1) is a bending dominated structure, 
built from three perpendicular crossed rings (designated as Circular lat- 
tice). In a previous study [26] on the cell geometry effect on in-plane en- 
ergy absorption behaviour of 2D honeycomb structures, it was shown 
that the octagonal honeycomb could provide the optimum energy ab- 
sorption performance after the hexagonal honeycomb. As hexagonal 
cell geometry was not easy to pack in three dimensions, and so an octag- 
onal cell was chosen for constructing 3D lattice in this work. Thus type 2 
was constructed from three crossed octagons, designated as Octagonal 
lattice. Its deformation mechanism under uniaxial compressive loading 
is a combination of bending and buckling of the cell edges. These two 
types of lattices (type 1 and 2) have been also investigated for tissue en- 
gineering scaffolds [6]. Some internal struts were added to strengthen 
the Octagonal lattice and form type 3 which was designated as 
Strengthen Octagonal lattice. It was shown that the vertical edges in 
the loading direction, which deform by buckling during the loading, 
cause instability and undulating behaviour of stress in the plateau 
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Table 2 


Material properties of HP PA12 determined from the standard tensile test (ASTM D638-14). 


(a) Physical and elastic properties 


Material Elastic modulus Poisson's ratio? Density” 0.2% offset yield strength 
HP PA12 1340 (MPa) 0.33 919 kg/m? 22.8 (MPa) 

(b) Post-yield true stress versus true plastic strain 

Yield stress (MPa) 22.8 27.6 31.3 35.6 40.9 476 51.0 51.7 
Plastic strain 0.000 0.007 0.014 0.025 0.045 0.085 0.129 0.150 


2 From FDM Nylon12 [37]. 
> Measured. 


region, which is not preferred for energy absorbing applications [26]. 
That's why the added struts in type 3 were inclined and not vertical. Lat- 
tice type 4 was built from tetrakaidekahedron cells, known as Kelvin 
cells, which is used for representing, idealizing and modeling of stochas- 
tic foams [27-29] due to its bending dominated deformation. Type 5 is 
constructed from rhombicuboctahedron (RO) cells and it has been in- 
vestigated for TE scaffold application [6,30,31]. Its deformation mecha- 
nism is stretching-dominated. Type 6 is constructed from simple cube 
cells which deform by buckling of the vertical struts under uniaxial 
compressive loading. 

All unit cell struts have a circular cross-section of diameter “d”. For 
unit cell of the Circular lattice type 1, “D” represents the diameter of 
the three crossed rings (refer to Table 1) and for all other unit cell 
types “I” represents the length of the equilateral members of polygons. 

Choosing this wide range of different unit cells allows the energy ab- 
sorption behaviour of lattices of different deformation mechanisms to 
be compared with each other. 


2.2. Fabrication of lattices using Multi Jet Fusion process 


In order to investigate the printability of the studied structures, a 
prototype of each lattice type was built using a new 3D printing technol- 
ogy called Multi Jet Fusion (MJF) in a brand of polyamide 12 plastic, 
namely HP PA12 available for this process. This recently developed tech- 
nique can build plastic parts with high quality, functionality and dimen- 
sional accuracy much faster than the currently available 3D printing 
technologies [32]. Another important advantage of this new process is 
the parts made by this technique are almost isotropic, and their me- 
chanical properties do not depend much on the building orientation. 

A schematic of the building process in MJF system is shown in Fig. 1. 
In this technique, the building process begins by applying a thin layer of 


powder material on the building platform. Then, and in one continuous 
pass, the fusing agent is applied to the material layer where the powder 
particles are meant to fuse together, and the detailing agent is applied to 
modify fusing and create fine detail and smooth surfaces. This pass also 
combines the printing with fusing energy, and reactions between 
agents and the material cause the material to selectively fuse together 
to form the part. 

Fig. 2 shows the printed prototypes of the investigated lattice struc- 
tures built by HP Jet Fusion 4200 3D printing machine in HP PA12 mate- 
rial. It was observed that the dimensions, appearance, quality and 
details of the parts were very good. 


3. Numerical simulation 
3.1. Material characterization 


To characterise the material properties of HP PA12, alongside with 
the lattices, four tensile standard test specimens of Type I with dimen- 
sions described in the standard test method for tensile properties of 
plastics, ASTM D638-14 [33], were built on HP Jet Fusion 4200 machine. 
The MTS Criterion Electromechanical Universal Test System (Model 43) 
with 50 kN load cell and an optical extensometer was used for the ten- 
sile test procedure at room temperature in accordance with ASTM 
D638-14. The resulting nominal stress-strain curves of the four samples 
under uniaxial tensile loading rate of 5 mm/min are shown in Fig. 3. The 
material behaviour of HP PA12 can be described as elastic, power hard- 
ening model [34]. It has relatively low yield stress but increases remark- 
ably in the plastic region. The Abaqus plasticity with isotropic hardening 
model was used to define the material's behaviour. The elastic proper- 
ties determined from the tensile tests are listed in Table 2a. The plastic 
nominal stress-strain data up to the necking obtained from the tensile 


Moving Plate 


Fig. 4. (a) FEA setup and boundary conditions of a lattice structure (b) compressing the lattice by moving the top rigid plate with a constant speed towards the fixed bottom rigid plate. 
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Fig. 5. Schematic compressive stress-strain behaviour of a cellular solid showing the three main regimes and the important parameters (reproduced from [37]). 


tests was converted to true stress-strain data using the Eqs. (1) and (2) 
given in [35] and is listed in Table 2b. These values are used in the FEA to 
define the elastic and post-yield behaviour of HP PA12 as required by 
the Abaqus material model [36]. 


Otrue = Onominal (1 + Enominal) (1) 


Etrue = In(1 + Enominal) (2) 


3.2. Finite element analysis of the lattice structures 


As the crushing mechanism of the lattices is a very complex process 
and involves all three types of nonlinearities, it cannot be simulated 
with implicit solvers. The three types of nonlinearities involved in the 
crushing of these structures include, firstly, the material nonlinearity, 
as the local stress in the cell struts exceeds the linear elastic limit of 
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the parent material and it undergoes plastic deformation (yields), sec- 
ondly, the geometrical nonlinearities due the large deformation of the 
structures, and the third type is the contact nonlinearity due to the in- 
teraction of struts during the crushing. 

Abaqus/Explicit package [38] was employed to simulate the lattices' 
compressive behaviour under quasi-static uniaxial loading. The material 
properties determined from the tensile test and described in Section 3.1 
were used as the material parameters input. Beam elements are the nat- 
ural choice for discretisation of such lattice structures. They are compu- 
tationally inexpensive and can be used to simulate models with many 
cells in a reasonable time [39]. Hence, the beam elements of type B31 
[38], a 2-node linear beam in space, were used to model the lattices. 
Through convergence studies, it was found that the element global 
size smaller than 1.125 mm could not improve the accuracy of results 
much while significantly increasing the computational time. Thus, all 
the lattice structures were discretised with the element size of 
1.125 mm. The compression simulation was performed using two 
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Fig. 6. Stress versus strain and energy absorption efficiency versus strain for (a) Circular and (c) Octagonal lattices, and their corresponding energy absorption diagrams for (b) Circular and 


(d) Octagonal lattices with plateau end indicated on the diagrams, all based on FEA. 
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(b) Strengthened 
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Fig. 7. Stress versus strain and energy absorption efficiency versus strain for (a) Strengthened Octagonal and (c) Kelvin lattices, and their corresponding energy absorption diagrams for 
(b) Strengthened Octagonal and (d) Kelvin lattices with plateau end indicated on the diagrams, all based on FEA. 


rigid plates at the top and bottom of the lattice structures. All degrees of 
freedom of the bottom plate were set to zero. The top plate was allowed 
to move just in the vertical direction and compressed the lattices when 
moving downwards at a constant speed towards the bottom plate, as il- 
lustrated in Fig. 4a and b. In order to prevent the penetration of 
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contacting surfaces during the crushing, a general contact algorithm 
(Explicit) was defined for the interaction of the lattice struts with them- 
selves and with the top and bottom plates. 

The compression simulation of each structure was performed up to 
the densification. The force and displacement data were obtained for 
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Fig. 8. Stress versus strain and energy absorption efficiency versus strain for (a) RO and (c) Cubic lattices, and their corresponding energy absorption diagrams for (b) RO and (d) Cubic 


lattices with plateau end indicated on the diagrams, all based on FEA. 
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Table 3 
Energy absorption and compressive properties of the six HP PA12 lattices derived from the 
FEA stress-strain diagrams. 


Type Lattice Plateau At the plateau end 


designation an Strain, Stress, Energy Efficiency, 
E Ecd o absorbed, E (%) 
(MPa) (MPa) W 
(J/cm?) 
1 Circular 0.79 0.68 0.99 0.51 51 
2 Octagonal 1.05 0.74 1.40 0.75 54 
3 Strengthened 0.85 0.67 0.96 0.55 57 
Octagonal 
4 Kelvin 0.89 0.71 1.10 0.59 53 
5 RO 0.67 0.70 0.96 0.46 47 
6 Cubic 0.48 0.85 1.46 0.41 27 


1000 evenly spaced time intervals during the loading cycle as suggested 
by ISO 13314:2011 [40]. 


3.3. Numerical stress-strain diagrams 


The compressive behaviour of cellular structures is commonly pre- 
sented as the compressive stress-strain diagram, in which values of 
compressive stress, O, are plotted as ordinates against corresponding 
values of compressive strain, €, as abscissas [41]. Here, the compressive 
stress, O, is the global stress and defined as the compressive load carried 
by the cellular solid specimen at any given moment divided by the orig- 
inal cross-sectional area of the specimen perpendicular to the loading 
direction expressed in force per unit area. The compressive strain, £, is 
defined as the dimensionless ratio of reduction in the gage length at 
any given moment to the gage length of the cellular solid specimen 
along the loading direction. 

Fig. 5 shows the schematic stress-strain behaviour of a typical cellu- 
lar material under compressive loading. When an open cellular solid or 
foam is compressed, initially the cell edges bend giving a linear elastic 
deformation. At a certain stress level, the cell edges begin collapsing 
by elastic buckling, plastic yielding or brittle fracture, depending upon 
the cell material properties. Subsequent collapsing of the cells continues 
at an approximately constant stress, termed plateau stress. When all cell 
edges collapse and contact each other, the densification regime begins, 
and the stress increases rapidly. At the end of this final stage, the densi- 
fied cellular solid behaves more like the dense parent material [3]. 

The compressive stress-strain diagrams of the investigated six lattice 
structures in this work under quasi-static displacement-controlled 
loading, as obtained from the FEA method described in Section 3.2 are 
presented in Figs. 6-8. The aforementioned three main regimes of the 
compressive behaviour of such structures, namely elastic, plateau and 
densification regimes, are apparent for each structure. Figs. 6-8 also 
show the efficiency versus strain and absorbed energy versus peak 
stress (energy absorption diagram) for all six structures that will be 
discussed in detail in later sections. 
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4. Energy absorption characterization of cellular structures 
4.1. Densification strain 


Plateau stress and densification strain are the most relevant proper- 
ties of a cellular structure to its energy absorption characteristic [34]. 
The plateau of the stress-strain curve of a cellular structure ideally ter- 
minates at the densification strain, €a. Theoretically, ¢g is equal to the 
porosity (1 — p*/p,) as it is the limiting strain at which all the pore 
space has been squeezed out. However, in reality, the cell edges jam to- 
gether at a smaller strain than this [1]. 

The optimum use of such structures is up to the plateau end (also re- 
ferred to as onset of densification), after which the stress increases sig- 
nificantly and become almost vertical at the densification strain [1], 
without much increase in energy absorption. 

The plateau end (onset of densification) strain and the densification 
strain relate to different stages of interaction between cell edges, and 
hence to different points on the compressive stress-strain curve (refer 
to Fig. 5). However many publications do not distinguish between 
these two critical strains [42]. It was shown that the method based on 
the energy absorption efficiency diagram (as discussed in Section 4.5) 
gives unique, precise and consistent results to determine the plateau 
end strain, &-q [42]. According to this method, the plateau end strain of 
each stress-strain curve was calculated as the strain that provides the 
highest energy absorption efficiency. Figs. 6-8 illustrate the plateau 
end of each stress-strain curve determined by this method and the pla- 
teau end strain, €¢q, for all lattices is listed in Table 3. Some other param- 
eters are listed in the table, which will be discussed in detail in the next 
sections. 


4.2. Plateau stress 


From Fig. 5 it can be observed that the amount of energy absorbed in 
the initial elastic regime is very small and most of the energy is absorbed 
in the plateau regime nearly at a constant stress level. Hence, the energy 
absorption capability and characteristic of any cellular structure mainly 
depend on the plateau stress magnitude and the plateau length. For a 
given relative density, the higher, the more stable (nearly constant) 
and the longer the stress plateau is, the better is the energy absorption 
performance. 

For a typical cellular solid, as shown in Fig. 5, the plateau regime 
starts from the crush (yield) strain, gy, representing the initiation of a 
new deformation mechanism of the cell edges, and ends at a critical 
strain, €ca, representing the end of stress plateau or onset of densifica- 
tion. The plateau stress is given by the equation [42], 


i o(e)de 
Ecd —Ey 


(3) 


Op = 


or simply it is the average of stress in the plateau region. For the studied 
lattice structures, the crush strain ey was determined as the first point 
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Fig. 9. Comparison of FEA efficiency parameter at the plateau end between the MJF HP PA12 (a) Octagonal and (b) Kelvin lattices. 


94 F.N. Habib et al. / Materials and Design 155 (2018) 86-98 


Stress, o 
Energy absorbed, 
Ww 


0% 20% 40% 60% 80% 100% 
Strain, € 


Stress, o 


100% 
80% 
60% 
40% 
20% 
0% 
0% 20% 40% 60% 80% 100% 
Strain, E 


Efficiency, E 


Fig. 10. Schematic (a) stress-strain, (b) energy absorption (energy absorbed vs peak stress) and (c) efficiency vs strain diagrams of an ideal energy absorber. 


on the stress-strain diagram at which an increase in strain occurs with- 
out an increase in stress if that occurred before the strain of 0.1. Other- 
wise, as in Circular and Kelvin lattices, the ey was taken as the strain of 
0.1 [41]. And, the efficiency diagram (Section 4.5) was used to deter- 
mine the plateau end strain, €.q. The calculated plateau stress for all 
studied lattice structures is listed in Table 3. Among all structures, the 
Octagonal lattice provides the highest plateau stress of 1.05 MPa. 


4.3. Energy absorption 


The energy absorbed per unit volume, W, by a cellular structure 
under compressive loading up to a strain £ is represented by the area 
under the compressive stress-strain curve up to the strain ¢ and given 
by [43]. 


W= 1 o(s)de (4) 


The cumulative area under each curve up to each 1000 strain inter- 
vals was calculated in Microsoft Excel using a numerical method, 
namely trapezoidal rule [44]. The energy absorption of each lattice at 
the end of plateau region, which is the optimum use of each structure 
[1] is listed in Table 3. From the table, it can be seen that the Octagonal 
lattice is capable of absorbing the highest energy of 0.75 J/cm?. 


4.4. Energy absorption diagram 


The energy absorption diagram [45] is a widely used engineering de- 
sign tool for characterizing the energy absorption behaviour of mate- 
rials and structures. It presents the relation between the two most 
important parameters for protective devices and packaging materials. 
It describes the amount of energy absorbed per unit volume, W, by a cel- 
lular material up to a strain € as a function of the corresponding peak 
produced stress, Op, on the protected object. From the diagram, it can 
be determined whether an energy absorber is capable of absorbing a 
particular amount of energy within the limited stress threshold to 
keep the protected object safe. 


Fig. 11. Mechanical compression test setup. 


The energy absorption diagram for the studied lattice structures is 
obtained by plotting the energy absorbed, W, against the peak produced 
stress, Op, up to each of 1000 strain intervals. The energy absorption di- 
agrams of the studied lattice structures are shown in Figs. 6-8. It can be 
observed that the plateau end, which was calculated as the point that 
provides the highest energy absorption efficiency, represents the shoul- 
der point on the energy absorption diagram precisely, which is consis- 
tent with the Maiti et al. [45] and Gibson and Ashby's [1] explanation. 


4.5. Energy absorption efficiency 


The efficiency parameter is the ratio of the energy absorbed by a ma- 
terial or structure up to a strain e to the energy absorbed by an ideal en- 
ergy absorber when both produce the same peak stress [46]. 

The efficiency parameter at the plateau end for the Octagonal and 
Kelvin lattices is illustrated in Fig. 9. It is equal to the area under each 
curve, W, up to the plateau end divided by the area under their corre- 
sponding ideal energy absorber (the green rectangular) which has a 
height of magnitude equal to the stress at the plateau end and a length 
of 1 (100% strain). 

The efficiency parameter, E, at a strain, £ is given by [46]: 


Op*1 
where op is the peak produced stress up to a given strain, €. 

In other words, the efficiency parameter, E, at any strain, £, is simply 
equal to the ratio of the energy absorbed to the peak produced stress op 
up to the strain €. 

Herein the ideal energy absorber needs to be defined. Such a struc- 
ture or material is the one which provides a rectangular stress-strain 
curve under compressive loading [34] as shown in Fig. 10a. An ideal 
structure does not deform initially (shows a rigid behaviour with no ini- 
tial elastic region) until the stress reaches the compressive strength of 
the energy absorber where the cells begin collapsing at a constant stress 
(plateau stress) up to 100% strain. It has the longest possible stress pla- 
teau, from the strain interval 0 to 1 without densification. In reality, 
there are foams made with a relative density as small as 0.001 [1], 
which theoretically means it can undergo up to 99.9% strain. The energy 
absorption diagram of an ideal energy absorber is shown in Fig. 10b. 

The schematic energy absorption diagram of an ideal energy ab- 
sorber is shown in Fig. 10b. It can be seen that there is no energy absorp- 
tion up to a certain stress (plateau stress), then the energy absorption 
goes up vertically at that constant stress. The efficiency of an ideal en- 
ergy absorber (Fig. 10c) has a linear relationship with strain. It starts 
from 0% at 0 strain and reaches 100% at 100% strain. 

The efficiency-strain diagrams for the lattice structures are shown 
on the stress-strain diagrams in Figs. 6-8 and the plateau end strain 
was determined as the strain which provides the highest energy ab- 
sorption efficiency as suggested by Li [42]. Among the six lattices, the 
Strengthened Octagonal and Octagonal lattices provide the highest en- 
ergy absorption efficiency of 57% and 54% as listed in Table 3. 

Some researchers have investigated the energy absorption capability 
of such materials using the so-called ideality parameter [26,43,46], 
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Fig. 12. Comparing the simulation (FEA) with the experimental (Exp) compressive stress-strain curves of HP PA12 lattices of (a) Octagonal and (b) Kelvin unit cells. 


which is the ratio between the energy absorbed by a cellular solid to the 
energy absorbed by an ideal energy absorber when both are com- 
pressed to the same strain, and produce the same peak stress [46]. How- 
ever, it was found that the ideality parameter cannot provide 
satisfactory results regarding the capability of cellular structures for ab- 
sorbing energy [43]. 


5. Experimental validation 
5.1. Compression tests of lattices 


To further investigate the compression response of lattice structures 
experimentally and validate the FEA compression results, two types of 
lattices (Octagonal and Kelvin) with the highest energy absorption pre- 
dicted by FEA method (Table 3) were selected for compression tests. 
Three specimens of each of these two lattices were tested under 
quasi-static uniaxial displacement-controlled loading. The same afore- 
mentioned universal testing machine (MTS Criterion Model 43) was 
used for the compression tests. The tests were performed at room tem- 
perature by placing a specimen on the lower fixed platform. The cross- 
head was moving downwards at a constant speed of 5 mm/min in 
accordance with the standard test method for compressive properties 
of rigid cellular plastics ASTM 1621-16 [47]. Fig. 11 shows the compres- 
sion test setup of an Octagonal lattice sample. The force-displacement 
data is recorded during the test on a connected computer. This data 
was converted to stress-strain data in the same way as explained in 
Section 3.3. The experimental stress-strain curves of three samples of 
both Octagonal and Kelvin lattices are compared to their FEA results in 
Fig, 12. 


(a) MJF sample 


It can be noticed from Fig. 12 that the FEA approach used is in good 
agreement with the experimental results. However, as reported in the 
previous studies [37,48] the FEA approach predicts a higher densifica- 
tion strain, much closer to the theoretical prediction [37] than the ex- 
perimental results, which causes some variance in the FEA and 
experimental results. The experimental energy absorption efficiency 
for the Octagonal lattice is found to be 48% versus 54% predicted by 
the FEA method, and for the Kelvin lattice, they were 49% vs 53% 
respectively. 


5.2. MJF process vs. FDM 


A prototype of each lattice also made by fused deposition modeling 
(FDM) process. Fortus 450mc FDM by Stratasys, which is one of the 
most powerful FDM systems currently available on the market, was 
used for printing the lattices. Careful attention was paid to choose ap- 
propriate process parameters to make sure the lattices would be printed 
properly. The lattices are constructed from cylindrical struts of small di- 
ameters ranging from 1.134 to 2.246 mm (refer to Table 1). To be able to 
print the cross-section of these cylindrical struts appropriately the 
smallest tip size (T12) available for Nylon12 (with the layer thickness 
of 0.1778 mm), and the smallest raster width of 0.254 mm were used. 
An Octagonal lattice sample made by this process (FDM) is compared 
to the same lattice made by MJF process in Fig. 13 and optical micro- 
graphs of a sample of each 3D printing process are shown in Fig. 14. It 
can be observed that the MJF part's appearance and quality are much 
better compared to the FDM sample. The time needed to print 4 lattices 
on FDM was 93 h and 16 min while for 12 lattices it took 3 h and 39 min. 
Which means the MJF required 0.013 fraction of the time required for 
the FDM process. The quasi-static compressive stress-strain curves of 


Fig. 13. An Octagonal lattice sample built by HP Jet Fusion 4200 (MJF) compared to a sample built by Fortus 450mc (FDM). 
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Fig. 14. Optical micrographs of (a) MJF sample and (b) FDM sample. 


MJF vs. FDM Octagonal lattice samples are compared in Fig. 15. It can be 
seen that the MJF sample provides much higher modulus, crush 
strength and energy absorption than the FDM sample. It was also ob- 
served that the FDM samples tear apart after the compression tests 
whereas the MJF samples showed no signs of rupture in the lattices 
struts after being fully compressed, and they even had good recovery 
behaviour, as shown in Fig. 16. 


6. Discussion of results 


From the comparison of MJF with the FDM lattices and their com- 
pression behaviour it was observed that the newly developed MIF pro- 
cess, could produce the lattice structures with much higher quality and 
functionality in 0.013 fraction of time needed for FDM process. 

The analysis of the FEA data, compression response and energy ab- 
sorption behaviour of the lattices showed that the lattices which deform 
by bending under the compression loading such as Circular, Octagonal, 
Strengthened Octagonal and Kelvin lattices generally provide good en- 
ergy absorption characteristic (high energy absorption and efficiency). 
Among the six studied structures, the Octagonal lattice is identified as 
the optimum structure with the closest performance to the ideal energy 
absorber. It provides the highest plateau stress of 1.05 MPa, the highest 
energy absorption per unit volume of 0.75 J/cm*, and the second highest 
efficiency of 54% as listed in Table 3. The high energy absorption perfor- 
mance of the 3D Octagonal lattice was found to be consistent with the 
results of 2D Octagonal honeycomb reported in [26]. The compressive 
response of the RO type lattice is stretching dominated and the Cubic 
lattice is buckling dominated, which makes them stiff and strong but 
unsuitable for energy absorption due to the significant softening of 
these structures after the yield point. These two lattices have the 
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Fig. 15. Quasi-static experimental (Exp) compressive stress-strain curves of MJF vs FDM 
Octagonal lattice samples. 


minimum energy absorption of 0.46 J/cm? and 0.41 J/cm? and minimum 
efficiency of 47% and 27% respectively. 

It was observed that when keeping the relative density constant, 
adding extra inclined struts to strengthen the Octagonal lattice of type 
2, caused an increase in efficiency from 54% to 57%. This is due to the 
more stable stress plateau of the Strengthened Octagonal compared to 
the Octagonal lattice. However, the energy absorption dropped signifi- 
cantly from 0.75 to 0.55 J/cm’, caused by a shorter plateau. 

The nonlinear FEA method used in this study provided predictions 
which were in good agreement with the experimental test results ob- 
tained from two types of lattices built by MJF. However, the FEA results 
predicted a slightly higher densification strain, and thus longer stress 
plateau, which resulted in overestimating energy absorption and effi- 
ciency. The methodology adopted in this study to determine the plateau 
end (onset of densification) by using the efficiency method, was pre- 
cisely consistent with the shoulder point on the energy absorption dia- 
gram as reported by Maiti et al. [45] and Gibson and Ashby [1]. The 
Octagonal and Kelvin lattices investigated experimentally in this study 
were found to provide higher energy absorption efficiency, 48% and 
49% respectively, when compared to the Octet-truss and all other cellu- 
lar materials reported in [25] including the Aluminum foam and 
honeycomb. 


7. Conclusions 


This investigation aimed to develop a high energy-absorbing lattice 
with the closest behaviour to the ideal energy absorber. The compres- 
sive stress-strain response and energy absorption characteristics of lat- 
tice structures from six different types of unit cells with identical 
relative density were evaluated and compared. A recently developed 
3D printing process, namely Multi Jet Fusion (MJF), was successfully 
used to print high-quality prototypes of the lattices in a brand of poly- 
amide (nylon) 12. Full-scale nonlinear FEA was used to investigate the 
compressive behaviour of lattices. The mechanical tests performed on 
two types of the lattices, Octagonal and Kelvin lattices, were in good 
agreement with the FEA predictions. Clear trends were observed on 
the influence of unit cell configuration on plateau stress, energy absorp- 
tion and efficiency. In terms of plateau stress, the lattices with bending- 
dominated deformation were found to be compliant and display a rela- 
tively stable plateau region resulting in higher energy absorption and ef- 
ficiency. On the other hand, those with stretching and buckling 
dominated deformation were found to be much stiffer and stronger, 
but soften after yielding, thus causing a lower energy absorption and 
efficiency. 

The results have identified that under quasi-static loading the Octag- 
onal lattice offers the optimum energy absorption characteristic among 
the six types of studied unit cell structures. It provides a high plateau 
stress, energy absorption and efficiency of 1.05 MPa, 0.75 J/cm? and 
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(a) MJF Octagonal sample 
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Fig. 16. Post-compression test samples built by MJF compared to an FDM sample. 


54% respectively. The results of this study reveal a strong relationship 
between the morphological and energy absorption properties of the lat- 
tice structures. The findings give a greater understanding on designing 
and 3D printing of lattice structures with the closest behaviour to the 
ideal energy absorber and can help to develop advanced lattice mate- 
rials that are superior to the currently used stochastic foams. This 
study focussed only on using the quasi-static displacement loading to 
examine the energy absorption behaviour of the studied structures. 
However, in impact applications, these structures will be exposed to 
various dynamic loading conditions depending on the application, 
which may cause a difference in their energy absorption behaviour. 
Hence, further studies are needed to investigate their performance 
under specified dynamic loading conditions. 
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